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Abstract. Evolution of galaxies is one of the most actual topics in astro- 
physics. Among the most important factors determining the evolution are two 
galactic components which are difficult or even impossible to detect optically: 
the gaseous disks and the dark matter halo. We use deep Hubble Space Tele- 
scope images to construct a two-component (bulge -I- disk) model for stellar 
matter distribution of galaxies. Properties of the galactic components are de- 
rived using a three-dimensional galaxy modeling software, which also estimates 
disk thickness and inclination angle. We add a gas disk and a dark matter halo 
and use hydrodynamical equations to calculate gas rotation and dispersion pro- 
files in the resultant gravitational potential. We compare the kinematic profiles 
with the Team Keck Redshift Survey observations. In this pilot study, two 
galaxies are analyzed deriving parameters for their stellar components; both 
galaxies are found to be disk-dominated. Using the kinematical model, the gas 
mass and stellar mass ratio in the disk are estimated. 
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1. INTRODUCTION 

Understanding the evolution of galactic structures is of uttermost importance 
in astrophysics. The knowledge of the processes determining the morphology, star 
formation history and intrinsic kinematics of galaxies give us the clues about the 
general cosmological framework: fundamental cosmological parameters, properties 
of dark matter, cosmic recycling and enrichment of baryonic matter, etc. However, 
several key aspects of galaxy evolution are poorly understood so far. 

Owing to extensive surveys, like the Sloan Digital Sky Survey (York et al. 
2000), the general structure of galaxies is quite well established (e.g., Simard et 
al. 2011; Lackner & Gunn 2012). For a direct tracing of evolutionary effects, the 
local galaxy sample has to be compared to galaxies at cosmologically significant 
distances. Unfortunately, only small samples of distant galaxies have been studied 
so far. From photometric decomposition of galaxies out to redshift z ~ 6 (e.g., 
Tamm & Tenjes 2006; Fathi et al. 2012), it has been shown that sizes of galaxy 
disks decrease with the redshift, in general accordance with simulations. 
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On the other hand, studies of the evolution of the intrinsic kinematics (and 

thus the gravitation potential and dynamical formation history) of galaxies do not 
extend much beyond z ~ 1. It has been shown that the evolution of the Tully- 
Fishcr relation is mild or missing over this period (Fernandez Lorenzo et al. 2010; 
Miller et al. 2011). 

In more detailed studies the photometric and kinematic data are used self- 
consistcntly to split the striicturc of galaxies into contributions by individual stel- 
lar, gaseous and dark matter components (Tenjes et al. 1994, 1998; Tempel & 
Tenjes 2006; Chemin et al . 2011; Jardel et al. 2011). However, for higher rcdshift 
objects the observational data rarely provide possibilities for such analyses (Tamm 
& Tenjes 2003, 2005). 

In the present study we have used both photometric and kinematic observations 
to study the structure of distant disk-dominated galaxies. We have applied a simple 
bulge + disk model on deep Hubble Space Telescope observations. Using the data 
for gas kinematics, a gas disk component was added to the model. Masses of the 
components were derived by solving the isotropic Jeans equations. 

Throughout this paper we assume a Friedmann- Robertson- Walker cosmological 
model with the total matter density fim = 0.27, dark energy density f^A = 0.73, 
and the Hubble constant Hq = 71 kms~^Mpc~^. 



2. DATA 



Photometry of galaxies is 
based on the deep Hubble Space 
Telescope ACS camera observa- 
tions in the GOODS project (Gi- 
avalisco et al. 2004). Modelhng 
is based on z filter observations, 
in which young stars and dust 
have the smallest effect. Initial 
data reduction, including dither- 
ing (Mutchler & Cox 2001) to 
0.03arcsec/pixel had been done 
by the GOODS team. Inten- 
sity counts were transformed into 
magnitudes using the formula 




1 arcsec 



1 arcsec 



Fig. 1. Photometric images of two galaxies 
(left TKRS 2568, right TKRS 9059) in z filter, 
with the position of the spectrograph slit over- 
plotted. 



m = -2.51ogio(^) + ^P, 



(1) 



where N is electron counts per second in CCD pixel and ZP is zero point of the 
used filter. According to Sirianni et al. (2005), the zero point for the z filter is 
24.86 mag. We took 4.51 as the value of solar luminosity. 

Figure 1 shows the z filter images of the studied galaxies, with the spectrograph 
slit positions and sizes indicated. Both galaxies have been selected to be disk- 
dominated objects with relatively good spectroscopic observations. 

We used spectroscopic observations from the Team Keck Redshift Survey (TKRS; 
Wirth et al. 2004) of galaxies in the GOODS field using the DEIMOS spectro- 
graph of the Keck II telescope. They have measured the [O II] 3727 A emission line 
along the slit roughly parallel to the major axes of galaxies. The slit width was 
1", which is comparable with the galaxy size. During spectroscopic observations 
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the seeing was 0.6" to 1". 

The original sample contains the rotation curves and velocity dispersion profiles 
of 380 galaxies extracted by Weiner et al. (2006) from the TKRS survey data. 
Since edge-on galaxies are difficult to use for photometric analysis and the face-on 
ones for kinematic analysis, we chose galaxies with intermediate inclination angles. 
Additional selection was based on the symmetry, spatial and velocity extent of 
the observed kinematics, required for reliable modeling. Approximately 10% of 
the galaxies matched these criteria. In this study we present the analysis of two 
galaxies (TKRS 2568 and TKRS 9059) which have good rotation and dispersion 
profiles and have also regular and disk-dominated morphology. The nomenclature 
is the same as used by Weiner et al. (2006). 



3. MODEL 

We used a three-dimensional model to describe the structure of each galaxy, briefly 
reviewed below. Details of the general model are given in Tempel et al. (2010, 
2011). 

3.1. PHOTOMETRIC MODEL 

The model galaxy is given as a superposition of its individual stellar components. 
In the present study we used two component model: bulge + disc. The light dis- 
tribution of both components are approximated by the Einasto law, 

l{a) = l{0) exp [- (a/fcao)^^^] , (2) 

where l{0) = hL/iA-Kqao) is the central density and L is the component luminosity; 
ao is the mean harmonic radius and a = sjR"^ + z^jq, where R and z are cylindrical 
coordinates and q is axial ratio of a component. The coefficients h and k are 
normalising parameters, dependent on N . 

These components are projected onto the plane of the sky using the following 
formula: 

L(X,F) = 2y^ riM^, (3) 

where = cos^ i + q"^ sin^ i is the apparent axial ratio, i is the inclination angle of 
a galaxy, A = \J X"^ + Y"^ /Q'^ is the major semi-axis of the equidensity ellipse of 
the projected light distribution. The summation is made over galaxy components. 

We estimated the correctness of the model fit, using the value, defined by 
the sum of the squared differences of the model and observational pictures. 

For minimizing we gave a physically justified initial guess for the parameters 
ourselves. In the next step, the parameter values were specified using the downhill 
simplex method of Nclder and Mead from the Numerical Recipe library. 

The photometric model gave us the stellar component (bulge and disc) param- 
eters (oo, q, N , L) and the inclination angle of a galaxy. These parameters were 
kept fixed during the following dynamical modelling. 

3.2. DYNAMICAL MODEL 

The dynamical model is based on the density distributions of the stellar compo- 
nents (bulge and disk), the gas disk and the dark matter halo. Density proffies of 
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the stellar components have the same distribution as the luminous components; 
a constant mass-to-light ratio was assumed, different for bulge and disk. The gas 
disk and dark matter parameters are derived during dynamical modeling. For 
describing the gas disk density distribution, we used Eq. ©. For the dark mat- 
ter distribution, we used the Einasto profile in another form, which is becoming 
increasingly popular for this purpose: 



PEinasto = Pc exp • 



(r/re)^/"-l \, (4) 



where n is in principle a free parameter. According to iV-body simulations, we 
take n = 6.0 (Merritt et al. 2006; Navarro et al. 2010). The term d„ is a function 
of n in a way that pc is the density at Tc defining a half-mass radius. The value of 
dn is 17.67 for n = 6 (Merritt et al. 2006). 

Dynamics of a galaxy was calculated from the potential and density distribu- 
tions using the Jeans equations. For simplicity, axisymmetry and stationarity of 
the galaxy were assumed for solving the Jeans equations. Since the observed dy- 
namics is based on gas motions, we could use the isotropic Jeans equations. The 
Jeans equations were thus used in the form 

where $ is the sum of potentials of galaxy components, V^p is rotational velocity, 
a is dispersion of velocities and p is gas disc density distribution; R and z are 
cylindrical coordinates. 

From these equations rotational velocity and velocity dispersions were derived 



V^iR,z) = -y 



dR ^dR 



(6) 



OO 



o^{R,z) = p~' j p{R,z')^^^^dz'. (7) 

To make the modeled data comparable with observations, we needed to find 
the line-of-sight projected velocity distributions. For constructing the line-of-sight 
distribution we applied a Gaussian velocity profile (based on V^p and a) in every 
point along the line of sight, and integrated them using the square of gas density 
as a weight. The projected velocity profile corresponding to each line 

of sight {X, Y) was found using the equation 

X 



Y , Vi?2 - X2 
sm I tan i 

where lgas{R,z) denotes the spatial gas density (given by Eq. ^) and f(R^z)iV) 
is the corresponding line-of-sight Gaussian velocity profile at a given point (i?, z) 
in the galaxy, calculated on the basis of Eqs. (O and ([7]). 
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Table 1. Parameters of galaxy models. 



ID 


z 


i 

[deg] 


Comp 


ao 
[kpc] 


1 


N 


itotal 


M 
[lO^Mo] 








Bulge 


0.9 


1.0 


5.1 


2.4 


7.0 


2568 


0.488 


49 


Stellar disc 


6.68 


0.2 


0.5 


22.3 


7.0 








Gas disc 


8.0 


0.1 


0.1 




1.0 








Bulge 


1.2 


1.0 


6.0 


3.7 


4.8 


9059 


0.253 


65 


Stellar disc 


5.95 


0.2 


0.7 


19.1 


7.0 








Gas disc 


10.0 


0.1 


0.2 




3.0 



3.3. COMPARING THE VELOCITY PROFILE WITH OBSERVATIONS 

The observational data points of velocity dispersion and rotation are given 
along the spectroscopic slit. To compare these data with our model, we had to 
take into account the width of the slit and seeing. To achieve this we firstly 
calculated the model output as a pixel map, where each pixel represents a line-of- 
sight integrated velocity distribution profile. 

The atmospheric effects were taken into account by convolution of the velocity 
profile map with a Gaussian seeing function. During the observations the seeing 
changed from 0.6" to 1" , therefore we used a Gaussian seeing kernel with a width 
of 0.8". Since the spectroscopic slit was of similar width (1"), the accuracy of the 
seeing kernel was not very important. 

Trying to take into account the slit width, we stmimed the velocity profiles 
along the slit in bins, centered to the given observational points. We used the 
square of spatial density of gas as a weight for the siimmation. Thus the final 
model velocity profile allows a direct comparison to the observational points. 

For fitting the model galaxy to the observed kinematics, we used the structure 
of the bulge and disk as derived during the photometric analysis. We used the 
maximal disk method, ascribing as much mass to tlic disk as allowed by the ob- 
served rotation curve. Subsequently, the masses of the bulge, gas disk and dark 
matter halo were derived to achieve the best fit to the observed kinematics. 

4. RESULTS 

Using the IRAF task ellipse, we found cUiptically averaged one-dimensional 
profiles of the model and the observed images using the same ellipse parameters 
for both cases. Figure 2 shows the model and the observed profiles in the z filter, 
together with the contributions by the bulge and disk components. The modeled 
galaxy parameters are given in Table 1. The results of the dynamical modeling 
arc presented in Figure 3, showing the observed kinematic profiles together with 
the modeled ones. 

Due to degeneracy, the found parameters are one of the possible parameter sets 
that fit the galaxy. For this the parameters of the dark matter halo arc responsible, 
because they originate from the inner parts of the galaxy, where the potential is 
dominated by the luminous components. 

TKRS2568 is a disk-dominated (bulge to total luminosity ratio B/T = 0.1) 
galaxy at redshift z = 0.488. The inclination angle is 49° , which makes the 
galaxy suitable for kinematic and photometric analysis. It has clumpy disk struc- 
ture, which could refer to a forming bulge (Conselice 2003). As can be seen from 
Figures 2 and 3, the four-component model (bulge + disk + gas + dark matter 
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Fig. 2. Results of photometric modelling for the z filter. 
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Fig. 3. Results of dynamical modelling for: TKRS 2568 and TKRS 9059. The x axis 

is from the slit center. 



halo) gives a good fit to the photometric and kinematic observations. Since the 
observed rotation curve does not reach the expected plateau supported by the 
gravitational potential of dark matter halo, the parameters for the dark matter 
component arc highly uncertain. 

TKRS 9059 is also disk-dominated galaxy (B/T = 0.14) at redshift 0.253. 
In the outer regions of the galaxy stellar streams or tilted spiral arms, maybe 
be re- lated to a recent merger or accretion event, are seen. They make the 
interpretation of the kinematic data less reliable: stationarity might not be fully 
achieved, which is also scon from higher dispersion errors. The inclination angle 
(67°) is determined from the inner, expected disk-dominated region. A good fit to 
the observed kinematics does not require the use of a dark matter halo component. 
Although not justified by the current cosmological paradigm, the data do not allow 
us to give any constraints to the possible dark matter halo. 

In general, resulting from the model constraints and the limited depth of obser- 
vations, the derived galactic components are somewhat degenerate and the pres- 
ence of the dark matter halos remain not confirmed. However, the modeling soft- 
ware enables to brake the degeneracy between the thickness and the inclination 
angle of stellar disk, and a corresponding study of larger samples would be impor- 
tant for detection of evolutionary effects in disk properties. 
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